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ABSTRACT 


Since the introduction of parameter plane techniques by Siljak 
in 1964, much effort has been devoted to the application of these 
techniques to various control system problems. 

This paper concerns the application of parameter plane techniques 
to the design of active filter systems. 

Constant bandwidth curves on the parameter plane are compared with 
frequency response curves to validate the constant bandwidth curves as 
accurate representations of active filter performance. This in turn 
implies that the constant bandwidth curves can be used in the design of 
active filter systems. 

Also presented is an example of the use of constant zeta and omega 
curves on the parameter plane to manipulate the pole locations of 
individual circuits to achieve a desired overall system configuration 


with the proper frequency scaling. 
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1. Introduction 

The use of operational amplifiers in active filters increases the 
flexibility with which signal filters can be designed. Primary reasons 
for this increased design flexibility are the elimination of the need 
for inductors and of the need for impedance matching. Active filter 
circuits can be used over a wide frequency range, oftentimes with the 
added benefits of less cost, smaller size and less weight. [2] These 
features suggest that active filters provide a suitable answer to a wide 
variety of filtering problems, such as, aerospace medicine data collec- 
tion and the analysis of vibrational fatigue in aircraft structures. 

This paper 1S a case study based upon an active filter parallel-T 
filter designed by Mr. William Kerwin of NASA's Ames Laboratory and is 
primarily concerned with establishing that parameter plane techniques can 
be effectively applied to design of active filter systems. Figure 1 is a 
circuit diagram of the filter and Figure 2 illustrates the observed 
experimental filter performance of that circuit. The experimental per- 
formance will serve as the fundamental reference for the evaluation of 
the parameter plane design techniques used in this paper. 

Two general types of parameter plane curves will be considered: 
(1) constant bandwidth curves and (2) constant zeta and omega curves. 

Hollister developed equations for displaying curves of constant 
bandwidth on the parameter plane and defined a constant bandwidth curve 


as: [3] 


A constant bandwidth curve for G(j@p) = M is a curve drawn 
upon the parameter plane which specifies the relation between 
the parameters necessary if the transfer function G(s), which 
is a function of the parameters, is to have magnitude M at 
the real frequency i). 
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Nutting developed computer programs to plot the constant zeta and 
constant omega curves on the parameter plane. B If these parameter 
plane curves are frequency normalized, then frequency scaling must be 
performed when cascading several circuits together to form a filter 
System. A technique for proper frequency scaling to allow the manipula- 
tion of system pole locations will be demonstrated. Proper conversion of 
parameter plane root locations into correct frequency scaled values for 
the individual circuits permits the individual circuit transfer functions 
to be multiplied together to obtain the system transfer function which 
represents the physical system. 

In establishing parameter plane techniques as applicable to active 
filter design, the following approach will be made. 

First, a mathematical model for the filter system will be developed, 
that is, the transfer function for the system will be derived. 

Then, a digital computer program for the CDC 1604 will be written 
to evaluate the frequency response of the mathematical model of the 
System. This frequency response will be compared with the observed 
experimental response. If the performance of the mathematical model 
conforms with the actual experimental response then the mathematical 
model will be considered a valid representation of the physical system. 
Furthermore, each of the individual circuits which together comprise the 
overall system will be considered to be accurately represented by their 
respective individually derived mathematical models, i.e., transfer 
functions. 

The frequency response of the transfer functions for the individual 
circuits will then be considered accurate enough to serve as standards 


for the evaluation of the constant bandwidth curves on the parameter 


10 


plane. Good correlation between the constant bandwidth curves and the 
corresponding frequency response curves then implies that both of the 
curves accurately describe circuit performance and that the parameter 


plane can, therefore, be used as a design tool. 


1a 


2. Mathematical Modeling of the System 

Evaluation of the effectiveness of parameter plane techniques 
could be accomplished in several ways. One method would be the Hard- 
ware approach. Use parameter plane techniques to determine the values 
of selected parameters which will produce a given performance, then 
actually build the system and adjust the chosen parameters to the speci- 
fied values. Observe the experimental performance and compare with the 
desired performance to substantiate the validity of the design technique. 

A second method is simulation. In this case a mathematical model 
of the system is developed and the performance of the system is evaluated 
on either a digital or analog computer. When facilities are available 
for either analog or digital simulation these methods have the advantage 
of rapid evaluation and excellent flexibility. 

Even if the first method of hardware analysis is applied, the devel- 
opment of the parameter plane curves requires an accurate mathematical 
model of the system. 

Therefore the first step in design using the parameter plane is to 
develop the mathematical model of the System. 

2.1. Derivation of Transfer Functions Through Signal Flow Graph Analysis 

The mathematical modeling of the filter system shown in Figure 1 will 
be based upon the analysis techniques of the signal flow graph introduced 
by S.J." Mason: [5] The system will be divided into the three circuits 
indicated in Figure 1. Transfer functions for the individual circuits 
will be derived using signal flow graph analysis. The overall system 
transfer function will be the product of the individual circuit transfer 
functions. This overall system transfer function will be evaluated by a 


computer program to be presented in Section 3. 


9.1.1. Signal Flow Graph Analysis ofCircuit.1. 
C 
1 





Cc out 


Figure 3. Circuit Diagram for Circuit 1 





Figure 4. Signal Flow Graph for Circuit 1 
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Non-touching loops taken two at a time, si 
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2.1.2 Signal-flow Graph Analysis of Active Parallel-T with Resistive 
Feedback. 


The circuit diagram shown below is an active parallel-T network 


and represents both circuit 2 and circuit 3. 





Figure 5. Circuit Diagram, Active Parallel-T Network 
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Figure 6. Signal Flow Graph, Active Parallel-T Network 
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The loops of the signal-flow graph, rie 
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Non-touching loops taken three at a time, Lye 
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multiply numerator and denominator by Pt, then factor (p + 1) from 
of 


numerator and denominator 


«{p2 + 1 
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Consideration of load impedance, Z, let 
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then define 
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Substitute into equation (2.2-1) 
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2.2. Frequency Sealing the Transfer Function 
Whenever several circuits are cascaded it becomes necessary to 


frequency scale the complete network to some reference frequency ae 
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prior to evaluating the overall frequency response. 
In the RC elliptical function filters under consideration the 


“notch frequency" is 
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Frequency scaling the transfer function: 


For example, consider filter section 2 
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Frequency scaling the transfer function. 


Similarly for filter section 3 
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Frequency scaling the transfer function. 
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3. Evaluation of the Mathematical Model 

Does the mathematical model developed in Section 2 accurately 
represent the physical system? Can those transfer functions be used 
in developing the parameter plane curves?. 

To answer these queStions the digital computer program, Program 
SYS RESP, written in FORTRAN 60, was developed to evaluate system fre- 
quency response when the system consists of several active circuits in 


cascade. 
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The system frequency response plotted by Program SYS RESP when 
the previously derived, frequency scaled circuit transfer functions are 
inserted as data is shown in Figure 7. Since it can be shown that 3.18 
ke corresponds to a normalized omega equal to unity, comparison of this 
ne system response with the experimental filter performance 
recorded in Figure 1 clearly indicates that the mathematical model is 
valid, that is, the transfer function is correct. In addition these 
results also indicate that the program will accurately simulate the per- 
formance of the physical system. 

Validity of the system mathematical model implies: 

(1) that each of the circuits which together comprise the system 
evaluated above can in turn be represented by its respective transfer 
Function. 

(2) that the performance of these individual circuit transfer 
functions can now be used in turn to evaluate the constant bandwidth 


curves on the parameter plane. 
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Figure 7. Digital Simulation of System Frequency Response 


4. Active Filter Constant Bandwidth Curves on the Parameter Plane 

The computer program included in Appendix A was written in FORTRAN 
63 by Jean Bow, Programmer, Computer Facility, Naval Postgraduate School 
and modified slightly by the author to display constant bandwidth curves 
on the parameter plane. The program is based upon the constant bandwidth 
equations developed by Hollister and utilizes the transfer function for 
the active parallel-T network derived in Section 2.1.2. 

The constant bandwitdh curves with amplitude contours in decibels 
are shown on the parameter plane plots of Figure 8-1 through Figure 8-18. 
As noted earlier each parameter plane plot of the constant bandwidth 
curves pertains to a fixed value of Q- 

To obtain a sketch of the frequency response of the circuit, select 
the values of the two parameters on the abscissa and ordinate. This 
establishes an operating point on the parameter plane. Enter each plot 
at the selected operating point to determine magnitude, M, at the fixed 
value of Ww, represented by that plot. The frequency response curve can 
now be constructed by plotting M as a function of G, 

Is this technique valid, that is, is it accurate enough to be con- 
sidered as a possible active filter design technique? 

The computer program included in Appendix B was developed specific- 
ally to answer these questions. Written in FORTRAN 60, Program 
FREQRESP evaluates the frequency response of the parallel-T active fil- 
Eel. 

Figure 9-1 is the graphical output for this program with SK, the 
ratio of the circuit element values, set at 2.0 and C, the normalized 
capacitance, equal to 0.20. The amplifier gain is varied from 1.6 to 


5.4 to produce the family of curves shown, i.e., each curve represents 
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The following information applies to Figures 8-1 
through 8-18: 


X Variable "Amplifier Gain, K" 

X Seale l inch = 0.5 

Y Variable "Normalized Capacitance, C" 
Y Scale finch = 0a 
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Figure 8-1. Active Filter Constant Bandwidth Curves on Parameter 
Plane k = 2.0 @ = 0.10 
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Plane k = 2.0 G@ = 0.20 
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Active Filter Constant Bandwidth Curves on Parameter Plane 


k = 2.0 
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Pigure 8-8. Active Filter Constant Bandwidth Curves on Parameter Plane 
k = 2.0 G) = 0.80 
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Figure 8-9. Active Filter Constant Bandwidth Curves on Parameter Plane 
k = 2.0 @ = 0.90 
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Pigure 8-10. Active Filter Constant Bandwidth Curves on Parameter Plane 
k=2.0 @ = 0.95 
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Figure 8-11. Active Filter Constant Bandwidth Curves on Parameter Plane 
k=2.0 @ = 1.05 
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Figure 9-1. Active Filter Frequency Response Curves 
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Figure 9-2. Active Filter Frequency Response Curves 
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a specific value of amplifier gain. At small values of w and again 
for G@ >)>1.0 we observe the expected result that output increases as 
gain increases. But the distinguishing feature of this family of curves 
is that near the "notch frequency" of 1.0 a resonance peak is observed at 
some intermediate value of amplifier gain. 

To evaluate the constant bandwidth curves, two operating points were 
selected. The first operating point corresponds to Circuit 2 in Figure 


1.» For this circuit 


SK = 2.0 
C = 0.436 
K = 2.762 


The frequency response curves for this circuit are given in Figure 9-2. 
A comparison of constant bandwidth against the frequency response values 
ls provided in Table I. 

The second operating point selected corresponds to Circuit 3 of 


Figure 1 for which 


SK = 2.0 
C = 0.284 
K = 2.074 


The frequency response curves for Circuit 3 are displayed in Figure 9-3 
and tabulated in Table II. 

Other operating points could well have been selected, but further 
repetition would do little more than reemphasize the conclusions that 
are arrived at by considering Tables I and II. For it is clearly appar- 
ent that the constant bandwidth curves on the parameter plane are excel- 
lent predictors of circuit performance. 

It is therefore valid to use the constant bandwidth curves for the 


design of such active filters. 


a2 


EVALUATION OF ‘'CONSTANT BANDWIDTH" CURVES 


TABLE I - Circuit 2 (C = 0.436, K = 2.762) 


RADTAN AMPLITUDE FROM AMPLITUDE FROM 
FREQUENCY CONSTANT BANDWIDTH CURVES FREQUENCY RESPONSE CURVES 

w ty My 
Oil 9.0 db 9.0 db 
OZ 9.3 9.4 

ORs JE Oe 10.0 

0.4 12 Lives 

O~g5 F3e0 iWSe6 

0.6 ESS 2020 

OE 7 lS8e0 19.0 

0.8 6.0 6.0 

0.9 -4,4 —4 4 

0.95 -12.1 =e 

1.05 -14.7 -14.0 

8, -9.8 -9.0 

iW 2e -5./ =) 0 

1.39 -3./ = 35.0 

1.40 -2.5 -2.5 

W360 =a -l.l 

1.80 -0.4 -0.3 
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EVALUATION OF "CONSTANT BANDWIDTH" CURVES 


TABLE II - Circuit 3 (C = 0.284, K = 2.074) 


RAD TAN AMPLITUDE FROM AMPLITUDE FROM 
FREQUENCY CONSTANT BANDWIDTH CURVES FREQUENCY RESPONSE CURVES 


W Ag Ax. 
Oe 6.2 db 6.2 sdb 
Ol a4 6.5 625 
Sls: Ga Oe a 
0.4 Grol tot 
0-5 Lae, i 
Oro Ti Pas 7a 
O27 620 6.0 
Diets) 1.8 1.4 
0.9 -6.1 =Oee 
0.95 =1357 -14.0 
ily Os =—la7 ~15.0 
Daetal sl -10.6 -10.6 
a2 -6.3 ~6.5 
iy 30 -4.3 -4,3 
1.40 -3.0 -3.1 
TOU -1.8 -1.8 
1.80 +O0.1 -1.0 
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5. Active Filter Constant Zeta and Omega Curves on the Parameter Plane. 

The computer programs developed by Nutting were the basis for the 
constant zeta and omega curves displayed on the parameter plane of Fig- 
ure 10. The transfer functions used for this parameter plane plot are 
derived from the active filter sections that have been designated circuit 
#2 and circuit #3. 

By careful selection of the proper values of amplifier gain, K, and 
the normalized input capacitance, C, it should be possible to adjust the 
root locations to any desired value of damping ratio, zeta, and radian 
frequency, omega, displayed on the parameter plane of Figure 10. 

5.1 Pole Locations and the Design of the Low Pass Filter 

The ideal low pass filter which is frequency response optimized 
has a gain vs. frequency characteristic shown in Figure ll. 

Although such filter performance is not physically realizeable, 
there are many filters that provide acceptable approximations to the 
ideal. Two mathematically designed filters, the Butterworth and the 
Chebyshev, will be discussed here. 

The Butterworth filter is said to be optimally flat at zero fre- 


quency, with amplitude related to frequency by 


where 


[a| is the amplitude 

Ww. is the cut-off frequency 

n 18 the order of the filter transfer function 
This relationship is plotted in Figure 12. 


The frequency response of the Chebyshev filter shown in Figure 13 


35 


Normalized Input Capacitance, C 


Figure 10. 





Amplifier Gain, K 


Constant Zeta and Omega Curves on the Parameter Plane 
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approximates the ideal frequency response with a "constant ripple" re- 
Sponse in the pass band, confining the error to some specified value 
throughout the pass band. It should be noted that the Chebyshev filter 
places equal emphasis on all frequencies in the pass band, whereas the 
Butterworth filter emphasizes the near zero frequency behavior. 

To achieve the frequency response shown in Figure 12, the poles of 
the Butterworth filter transfer function displayed on the s-plane are 
uniformly spaced on a semicircle of radius @ as shown in Figure 14. 

In a similar fashion the frequency response of the Chebyshev 
filter is accomplished by placing the poles of the filter transfer 
function on an ellipse on the s-plane. It has been shown that the proper 
positioning of the poles can be related to the poles of a Butterworth 
Filter of the same order by reducing the real component of the respective 
roots by the factor tanh & [6], without modifying the imaginary component. 
Thus 


om 


Chebyshev 7 (DRutterworth) tanh #) 


where 


The 


i) SS =) |] 
as Sinh co 


n is the order of the filter 
€ is related to the amplitude of the ripple 


ee the minimum amplitude in the 


V1 + & 


pass band as shown in Figure 13. 
5.2. Manipulation of Circuit Poles into Chebyshev Configuration 
The poles of the active filter system shown in Figure 1 will now 
be modified utilizing parameter plane techniques to form a Chebyshev 


configuration of poles. 
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Figure 11. Ideal Low Pass Filter Frequency Response 
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Figure 12. Butterworth Filter Frequency Response 





Figure 13. Chebyshev Filter Frequency Response 
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Figure 14. Relative Position of Poles for Butterworth 
and Chebyshev Filters (6th Order) 
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The purpose of this exercise is two-fold: 

(1) To illustrate the methods required for proper pole relocation 
which must include an adjustment of frequency scaling factors introduced 
through the normalization of the notch frequencies that was one of the 
steps in producing the normalized parameter plane of Figure 10. 

(2) To demonstrate that the constant zeta and omega curves on the 
parameter plane are an effective means of obtaining a desired performance 
from an active filter system. 

For the purposes of this discussion the roots associated with the 
low pass filter circuit #1 will be assumed to be fixed and the positions 
of the remaining poles will be manipulated to obtain the Chebyshev con- 
figuration. Another assumption is that the notch or "zero output" fre- 
quency is high so that the effects of the transfer function zero will be 
ignored, that is, only the poles will be considered. The two assumptions 
are made for the convenience of clarifying the following presentation. 

No serious limitations upon the design of active filter systems are 
introduced by these assumptions. 

Figure 15 is an s-plane plot of the pole configuration of the active 
filter system. The nuinneee in parentheses identify the circuit number 
which is associated with that particular pole pair. Thus the poles 
nearest the negative real axis designated (1) are associated with cir- 
cuit #1 and form the basis for the placement of the remaining two sets 
of poles. 

Since the imaginary component of this pole pair is equal to the 
imaginary component of the corresponding Butterworth filter, the poles 
designated (1) establish the radius, Coe of the reference Butterworth 


semicircle. 
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Figure 15. Active Filter System Pole Configuration 


61 


QJ). sine165° =80 23638 


@ (1) 
pg) reno 
sin 165° 


Similarly the real component of the poles designated (1) establish 


the magnitude of the ripple by specifying the value of tanh 


tanhg - OT Chebyshev - o~ (1) 
O putterwrth wd, cos 165° 
tanh B - -0.5/776 = 0.4254 


(-0.9659) (1.406) 


Now that Q | and tanhA have been eStablished, the remaining pole 


positions can be determined. 


G) (3) = a. Sin 135° = 0.994 


1.358 


Q(2) = QJ). sin 105° 


= = O 
CO) aA nvanen - tanh Op utterworth a tanhg a, coors 


= -0.578 
pond O = = 
Se Che paetiev O. 4254 a cos 135 Oe423 
— . a O = _ o 
oo) 0.425 a) cos 105 0.155 


The values of the damping ratio for the respective roots, zeta, 
will be used on the parameter plane of Figure 10 because of their qual- 
ity of maintaining the same value completely independent of any frequency 
Scaling or frequency normalization that must now be performed in order 
to take advantage of the information contained on the parameter plane 
plot of Figure 10 which, it must be recalled, has been normalized to a 


"notch" frequency equal to 1.0. 
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To calculate the values for zeta the following relationship must be 


applied: 
a 
Vo-? ss DG 
Thus: 
#57 8 


io rey 
Cheb \ /0.334 + 0.132 


0.423 


((3) =o OEE ee 
Sass \/0.179 + 0.988 


OFM SYS 


(2) SS 
Gaels 
¢ = \ /0.024 + 1.843 


In summary, the following values will produce a Chebyshev config- 


0.847 


0.392 


Oy ts 


ULrat Lom: 
Cireuit #1 
(4) = 0.847 
@ (1) = 0.364 


oO (1) = -0.578 


€(3) = 0,392 
(3) = 0.994 


o7 (3) = -0.423 


A= 0.113 
GJ (2) = 1.358 
O (2) = -0.155 
Since the parameter plane of Figure 10 is normalized to a "notch" 
Frequency of 1.0, the calculated frequencies listed above must be divided 


by the transfer function scaling factors, a, discussed in Section 2.2. 
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The scale factors for circuit 43 and circuit 42 were 1.234 and 1.566. 


Therefore the parameter plane frequencies, @ , are: 
p 


_ 0.9941 
1.234 


0.8056 


od 63) 


Ue Seg 
2) Gia. = 08867 
a 1.566 : 


To obtain these parameter plane values of omega and zeta, enter 
Figure 10 to determine the required values of normalized input capaci- 
tance, C, and amplifier gain, K. The following parameter values will 
result: 


Cireult #73 


KC3) = 1.605 a(3) = 1.234 
C(3) = 0.100 

Girecult #2 
K(2) = 2.04 a(2) = 1.566 
C(2) = 0.107 


Refer to Appendix C for algebraic verification that a Chebyshev 
configuration has been formed by the parameter values selected above. 
5.3. Evaluation of Response of System with Pole Configuration Derived 

Through the Parameter Plane. 

What is the effectiveness of the parameter plane in realization 
of the "desired performance", which in this case was constant ripple in 
the pass band? Program SYS RESP, presented in Section 3 was used to 
conduct this evaluation. 

The overall system transfer function required as input data to 
Program SYS RESP demands proper frequency scaling of the individual 


circuit transfer functions according to the formula derived in Section 2: 
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K(p* + a2) 


iL. (poe 
7 (2Gelpe + a (3.0 © Lok) +30 p +a 


2 


Circuit #3 
Z 
. WeGsiGpe + 1.522 
1, 30p- @ telepe+ 1.522 
Circus #2 


2.04(p? + 2.449) 


T 9 (P) = 
2 1.321p + O0.408p + 2.449 


The fixed pole associated with circuit #1 retains its original 


elrcuit transfer function: 


-9316 


Th? SSeS 
p* + 1.1553p + .4658 


The results of the Program SYS RESP evaluation of the active filter 
system designed through the application of constant zeta and constant 
omega curves on the parameter plane are shown in Figure 16. These re-~ 
sults indicate that the design objective of obtaining constant ripple 


performance has been attained. 
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Figure 16. Digital Simulation of Chebyshev System Frequency Response 
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6. Conclusions 

It has now been established that the constant bandwidth curves on 
the parameter plane and the constant zeta and omega curves on the para- 
meter plane provide valid information which may be used in the design 
of active filter systems. 

The constant bandwidth curves can provide a very accurate indica- 
tion of filter performance which is directly related to parameter values. 
In addition the constant bandwidth curves provide the filter designer 
with a good qualitative indication of parameter values which should be 
avoided because of extreme sensitivity of filter performance to varia- 
tions in those specific parameter values. Such regions of increased 
sensitivity are easily recognizeable by the increased density of the 
contour lines on the constant bandwidth plots. 

When desired filter performance can be realized with a known pole 
configuration, the constant zeta and omega curves provide a very effec- 
tive means of achieving that desired performance through the proper 
specification of the various parameter values. These curves can also 
provide the user with information pertinent to the tuning of the filter 
when performance is affected by tolerances of the circuit elements or 
the aging of components. 

It should be noted, however, that the frequency normalized para- 
meter plane can yield meaningful information only with proper frequency 
scaling. Considerable care must be exercised to correctly frequency 
scale the transfer functions when several circuits and/or systems are 
cascaded to achieve some desired result. 

Subject to these conditions parameter plane techniques are useful 


tools in the design of active filter systems. The information contained 


6/7 


on various parameter plane presentations should lead to much greater 
insight into the problem of correlating parameter values, pole locations 
and filter performance. Used in conjunction with the tremendous flexi- 
bility and freedom that prevails in the relocation of active filter poles, 
these parameter plane techniques should be substantial assets in further 


development of filter design theory. 
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APPENDIX B 
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OXVIVON OC) ORO) Cy CNR 1 GO OC: oe 


FUNCTION HAVING A NUMERATC2? POLYNOMIAL OF THE FOR” 
K x (PK¥2 + 1) 
AND A DEVOAINATOR POLYNOMIAL GF THE FORM 


Rete +) ) eee 


(Come ye 02 —Kh)/SK4IUSK4+1)*C)#P 4 7 


WHERE 
SK IS THE RATIO OF ELEMENT vALUES 
C 1S THE NORMALIZED CAPACITANCE 
EACH PLOT WILL sBE FOR" mipG@iieeeileE OF eGK AND C, BUT 
A FAMILY OF CURVE? WILL RE PLATTED TO REPRESENT VARIOUS 


VAISeo 


DIwkNSTUN 
LBREAL( 590)» SIMAG( 5 GOT iro AREAL ( SU Gala WiAG(500)> 


OF AMPLIFIER 


GAIN Ke 
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~—20B(50U)s AMAGN(5U0) 
READ 2oG, (ITITLEC]) > 1l=iaay 
READ 20s (ITITLECI)s I> 7eetee) 
EAtI) = 41.6 
LA(2) = 4HK1.2 
ON a I 2 eC 
LA(4) = 4HK2.2 
AS = 4K 2eg4 
LA(6) = 4SHK2.6 
LA(7) = 4HKZ2.8 
ewes = 4k 
LA(9) = 4HK3.2 
LA(1U) = 4HK3e4 
EACl i) SHKBe6 
LA(12) = SHK3.-8 
LA(13) = 4HK4.C 
LA(14) = 4HK4.2 
LA(15) = 4HK4.4 
LA(16) = 4HK4-e6 
SK = 620 
C = Ueb60 
MAXK = 16 
GAIN = 1-4 
w(l) = Q.2000 
NUMPTS = 200 
DO 1GQ0 I = ls MAXK 
GAIN = GAIN + O42 
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22 GSN TINUE 
ee (=>) 33 9557550 
33 MODE = 1] 
SOTO -55 
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37 MODE = 3 
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200 FORMAT (6A8) 
END 
END 
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APPENDIX C 

Algebraic Verification of Chebyshev Pole Configuration 

The parameter plane techniques of Section 5 should yield parameter 
values which will result in the desired Chebyshev pole configuration. 
An algebraic evaluation of the characteristic equations for circuits l, 
2, and 3 will now be conducted to confirm that such a Chebyshev pole 
configuration has, in fact, been attained. 

The parameter plane produced the following frequency scaled transfer 


functions for circuits 1, 2 and 3 respectively: 


0.9316 
A). ae ie 
p~ + 1.1553p + 0.4658 
2.04(p2 + 2.449) 
Too (PD = Ce ees aa 
1.321p* + 0.408 + 2.449 
1.605(p* + 1.522) 
T .(p) = => ar 
o3 1.30p2 + 1.100p + 1.522 


The characteristic equations (denominator polynomials) for these 


transfer functions are of the form 
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i} 
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ap + Bp er © 


since 
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and the imaginary 


a) 


Cireurt: | 


@® (1) 
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(1) 


CU (1) 


Gireurt 2 


(J) (2) 
n 


((2) 


C0) (2) 
Circure 3 


G (3) 
n 


C(3) 


(3) 


7 Bw 


n 
2G 


component of the pole becomes 


w@ \f1 -¢° 


0.4658 = 0.6825 


1.155) (0.6825 
e825) ee 
2(0.4658) 


(0.6825) (0.5331) = 0.364 


2.449 
- eS on LienG2 
15321 


_ _(0.408) (1. 362) 
2 (2.449) 


= Oebe2 


Yr. 354 


reso2) (079936) = 


ge 1.082 
Pao) 


(1.100) (1.082) 
= Ee Tae =. 6 SURE 
DA 22)) “ 


(1.082) (0.9204) = 0.996 


Comparison of these 
desired Chebyshev values 
of the constant zeta and 


Sign task of obtaining a 


calculated values of zeta and omega with the 
of Section 5.2 substantiates the effectiveness 
omega curves on the parameter plane for the de- 


desired active filter system pole configuration. 


80 


10. 


INITIAL DISTRIBUTION LIST 


Defense Documentation Center 
Gameron station 
Alexandria, Virginia 22314 


Library 
Naval Postgraduate School 
Monterey, California 93940 


Dr. George Thaler 

Dept. of Electrical Engineering 
Naval Postgraduate School 
Monterey, California 93940 


Die ook. barker 

Dept. of Electrical Engineering 
Naval Postgraduate School 
Monterey, California 93940 


Dr. Dee Dee oak 
University of Santa Clara 
Sant a eUaraweGarlitornia 95053 


Mr. William Kerwin 
NASA Ames Aeronautical Laboratory 
Moffett Field, California 94035 


LT David Freese, USCG 
Eleventh Coast Guard District Hdqtrs. 
Long Beach, California 90802 


LCDR G. R. Nakagawa, USN 

Attack Squadron 165 

c/o Fleet Post Office 

San Francisco, California 96601 


Dr. L. Eisenberg 

Newark College of Engineering 
323 High Street 

Newark, N.J. 07102 


Do oR, Low Ut 

Reader in Engineering Production 
Welsh College of Advanced Technology 
Cathays Park 

Cardiff, Wales 


oil 


No. Copies 


20 





TAB ba ‘ay 
| ii} if 3 (- 
, Oy 











"i 
if! 
| 


| f 





UNCLAS SIF IED 
Security Claaaification 


DOCUMENT CONTROL DATA - R&D 





(Security cleesification of title, body of abetract and indexing annotation muet be entered when the overall report ie claeeified) 
11. ORIGINATING ACTIVITY (Corporete author) Za. REPORT SKCURITY CLASSIFICATION 


Naval Postgraduate School UNCLAS SIF IED 


—— oe eT 


PARAMETER PLANE TECHNIQUES FOR ACTIVE FILTER SYSTEMS 







3. REPORT TITLE 





4. DESCRIPTIVE NOTES (Type of report and inclueive datee) 
Thesis, M.S., December 1966 


8. AUTHOR(S) (Laet name, fret neme, initiel) 


NAKAGAWA, Gordon R. 





8. REPORT DATE Ja. TOTAL NO. OF PAGES 75. NO. OF REFS 
December 1966 82 6 


Ba. CONTRACT OR GRANT NO. 9a. ORIGINATOR'S REPORT NUMBER(S) 
b. PROJECT NO. 


9 b. OTHER Ryroe NO(S) (Any other numbere that may be aeeigned 
fe repo 


10. AVAILABSILITY/LIMITATION NOTICES 


11. SUPPL EMENTARY NOTES 12. SPONSORING MILITARY ACTIVITY 


13. ABSTRACT 


Since the introduction of parameter plane techniques by Siljak in 
1964, much effort has been devoted to the application of these techniques 
to various control system problems. 


This paper concerns the application of parameter plane techniques to 
the design of active filter systems. 


Constant bandwidth curves on the parameter plane are compared with 
frequency response curves to validate the constant bandwidth curves as 
accurate representations of active filter performance. This in turn 
implies that the constant bandwidth curves can be used in the design of 
active filter systems. 


Also presented is an example of the use of constant zeta and omega 
curves on the parameter plane to manipulate the pole locations of 
individual circuits to achieve a desired overall system configuration 
with the proper frequency scaling. 





DD .52%%. 1473 UNCLASSIFIED 


Security Classification 


UNCLASSIFIED 
Securi assification 


KEY WORDS 


active filters 
active networks 
Butterworth filter 
Chebyshev filter 
filters 

frequency response 
networks 

parameter plane 


Proce] wr [nove | wr [nove | wr 


INSTRUCTIONS 


1. ORIGINATING ACTIVITY: Enter the name and addresa 
of the contractor, subcontractor, grantee, Department of De 
fanse activity or other organization (corporate euthor) iasuing 
the report. 


2e. REPORT SECURITY CLASSIFICATION: Enter the over 
all aecurity cleseification of the report. Indicate whether 
‘*Restricted Date’’ is included. Marking is to be in accord 
ance with appropriate security regulations. 


25. GROUP: Automatic downgreding is specified in DoD Di- 
rective 5200.10 and Armed Forces Industrial Manual. Enter 
the group number.- Also, when applicable, show that optional 


markigga have been used for Group 3.and Group 4‘as author- . | 
ized. 


& as 3 ae rae 
3. REPORT TITLE: Enter the complete report title in all 
cepital lettere. Titles in all ceses should be unclesailfied. 
If a meaningful title cannot be selected without classifice 
tion, show title clessification in all capitals in parenthesis 
immediately following the title. 


4 DESCRIPTIVE NOTES. If appropriete, enter the type of 
report, eg, interim, progress, aummary, annual, or final. 
Give the inclusive detes when a specific reporting period is 
covered. 


5S. AUTHOR(8® Enter the name(s) of author(s) as shown on 
or in the report. Eater leat name, first name, middle initial. 
If xilitary, show rank and branch of aervice. The name of 
the principal author is an absolute minimum requirement. 


6. REPORT DATE: Enter the date of the report as day, 
month, year, or moath, yean If more than one date appears 
on the report, use date of publication. 


7a. TOTAL NUMBER OF PAGES: The total page count 
should follow normal pagination procedurea, i.e, enter the 
number of pages containing information. 


76. NUMBER OF REFERENCES Enter the total number of 
referencea cited in the report. 


Be. CONTRACT OR GRANT NUMBER: If appropriate, enter 
the appliceble sumber of the contract or grant under which 
the report was written. 


85, &, & 66. PROJECT NUMBER: Enter the appropriate 
military depertment ideatificetion, such as project number, 
subproject auaber, system numbera, task number, etc. 


9a. ORIGINATOR‘'S REPORT NUMBER(S): Enter the offi- 
cial report aumber by which the document will be identified 
and controlled by the originating activity. Thia number muat 
be unique te this report. 

9b. OTHER REPORT NUMBER(S): If the report has been 
assigned eag ether report sumbera (aither by the originator 
or by the epeaecoe), alee enter this number(s). 


10, AVAILABILITY/LIMITATION NOTICES: Eater any lice 


itetions en fwther Giseemination of the report, other than these 


DD «18%. 1473 (BACK) 


& 


imnoane by security classification, using standard statements 
auch as: 


(1) ‘*Qualified requesters may obtain copies of this 
report from DDC.’‘ 


(2) “Foreign announcement and dissemination of this 
report by DDC is not euthorized. ‘‘ 


(3) ‘‘U. S. Government agencies may obtein copias of 
this report directly from DDC. Other qualified DDC 
users shall request through 


9 
° 


- (4) *U. S. military agencies may obtain copies of this 
.teport directly from DDC Other qualifjed users 
shall requeat through 7 
9 "CO wy 


eT take j 





o¢ 
° 





(S) ‘All distribution of this report is controlled Qual- 
ified DDC users shall request through 


e¢ 
d 


If the report has been furnished to the Office of Technical 
Servicea, Department of Commerce, for sale to the public, indi- 
cate thia fact and enter the price, if known 


11, SUPPLEMENTARY NOTES: Uae for additional explane- 
tory notes. 


12, SPONSORING MILITARY ACTIVITY: Enter the name of 
the departmental project office or laboratory sponsoring (pay 
ing for) the research and development. Includa address. 


13. ABSTRACT: Enter en abstract giving a brief and factual 
summary of the document indicative of the report, even though 

it may also appear elsewhere in the body of the technical re- __ 
port. If additional space is required, a continuation sheet shal! 
be ettached. 


It is highly desirable that the abstract of classified reporta 
be unclassified. Each paragraph of the abstract shall end with 
an indication of the military security classification of the in- 
formation in the paragraph, represented es (TS), (8), (C). or (UV). 


There is no limitation on the length of the abetrect. How- 
ever, the suggeated length is from 150 to 225 words. 


14. KEY WORDS: Key words are technically meaningful terms 
or short phrases thet characterize a report and may be used as 
index entries for cataloging the report. Key words must be 
selected so thst no security classification is required. Identi- 
fiers, auch as equipment model designation, trade name, military 
project code name, geographic location, may be used as key 





words but will be followed by an indicetion of technical con- 


text. The assignmant of links, reles, and weights is optional. 


UNCLASS IFIED 
Security Classification 


‘4 














nd Se 














